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Abstract 

We review and compare the reconstruction methods of the inclusive deep inelastic 
scattering variables used at HERA. We introduce a new prescription, the Sigma (E) 
method, which allows to measure the structure function of the proton F2{x, Q^) in a large 
kinematic domain, and in particular in the low x-low region, with small systematic 
errors and small radiative corrections. A detailed comparison between the S method and 
the other methods is shown. Extensions of the S method are presented. The effect of QED 
radiation on the kinematic reconstruction and on the structure function measurement is 
discussed. 



1 Introduction 



The measurement of the structure functions of the nucleon requires a precise reconstruction of 
the deep inelastic scattering (DIS) kinematics. With the advent of HERA, the first electron- 
proton collider ever built, this reconstruction needs not anymore to rely on the scattered 
lepton only, since the most important part of the hadronic system is visible in the An detectors 
HI and ZEUS. This redundancy allows for an experimental control of the systematic errors 
on the structure function measurement if it is based on several independent methods to 
determine the usual DIS kinematic variables x,y and Q^: 

= -(k - k'f = -q^ Q^ = xys (2) 

with s being the center of mass energy squared, P, k the 4-vcctors of the incident proton 
and lepton and k' the scattered lepton one. In this report we briefly review in section 2 
the methods used at HERA in the first year of operation. We introduce in section 3 a new 
way to compute the kinematics, the S method, which allows a precise F2 measurement in 
the complete kinematical region accessible at HERA including the low x-low region. We 
then make two rather detailed comparisons of the E, the mixed, the double-angle and the 
electron methods at low x and at high (sect.4). Extensions of the S method, including 
a prescription for kinematic fitting are given in section 5. We examine the effect of QED 
radiation on the F2 measurement (sect. 6). Two methods independent of initial state QED 
radiation are defined and compared in section 7. The main conclusions are summarized in 
section 8. 



2 Kinematic Reconstruction Methods at HERA 

In the naive quark-parton model, the lepton scatters clastically against a quark of the proton, 
and the two body final state is completely constrained using two variables, if we know the 
initial energies labeled and EP for the electron and proton. Denoting E,9 and Fq,jq the 
energy and angle of the electron and of the struck quark (which gives rise to the "current 
jet"), we can derive y and by combining any 2 of them. For a scattered quark Fq and 7^ 
are well defined observables, but what is actually observed after a collision is a set of particles 
making energy deposition in the detector. After identification of the scattered electron, we 
can define 2 independent quantities: S, defined as the sum of the scalar quantities Eh—Pz,h of 
each particle belonging to the hadronic final state, and T as its total transverse momentum. 

j: = Y.iE^-Pz,n) T= l{Y:P.,hr + i^Py,hr (3) 

h \ h h 

Eh,Px,h,Py,h,Pz,h are the four-momentum vector components of each hadronic final state 
particle. Prom energy-momentum conservation we can then define an inclusive angle 7 but 
also an inclusive energy F of the hadronic final state which would be the same as the angle 



1 



and energy of the scattered quark in the naive quark-parton model. The formulae which 
allow to go from S, T to F, 7 are: 



tan ^ = - F= (4) 

S is by construction minimally affected by the losses in the forward direction]^ (= 2Ei sin^ ^ 
for every lossed particle i) due to the beam pipe hole in which the target jet and the initial 
state gluon radiation tend to disappear. T covers the other spatial dimensions, but it is more 
sensitive to forward losses (= EisinOi), thus to improve kinematic reconstruction we should 
try to avoid using it. We can use 7 instead, which carries the T information and is better 
measured since in the ratio S/T the energy uncertainties cancel and the effect of the losses 
in the forward beam pipe is diminished. In conclusion the optimal four "detector oriented" 
variables to characterize deep inelastic scattering at HERA are [E,9,T,,^]. 

If we want to determine x and from 2 variables out of these 4, we are left with 6 com- 
binations of which only 3 have been shown to be usable for structure function measurement 
(the formulae are given in the appendix): 

• the electron-only (e) method which uses E and 9 is experimentally simple and very 
precise in the high y (> 0.2) region. It nevertheless degrades rapidly with decreasing 
y [01 due to the 1/y term in the error propagation^ of ye'- 

6y_l-y 6E 5^ 
y - y ^ E ®tane/2> 

• the hadrons-only (h) method |^] provides a rather precise measurement of y at low and 
medium y {y < 0.2), which degrades at high y. It gives a rather poor measurement of 

due to the loss of hadrons in the beam pipe but it is the only possible method for 
"charged current" events. 

• the double-angle (DA) method |3|, |^ uses the angle of the electron 9 and the inclusive 
angle of the hadronic final state 7. Assuming an homogeneous energy measurement 
over the full solid angle, we can deduce from eq.4 that the DA method is independent 
of the absolute energy calibration. This method is precise at high (Q^ > lOOGeF^) 
but has large resolutions at low x-low , as will also be illustrated in section 4. 

At HERA a fourth method has been used which combines y measured from the hadrons and 
from the electron, which is the so-called mixed (m) method |^: 

with yh = ^ (6) 



syh 2E- 

The mixed method allows to extend towards lower y the F2 measurement done with the 
DA method |^] or with the e method 0. Note that in the mixed method more than 2 

^The positive z axis is defined at HERA as the incident proton beam direction. The electron polar angle 
is thus large (~ 170°) for small scattering angle, i.e. small . 



We define A®B = ./WTB^ 
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independent variables are used, i.e. E, 6 and S, justly avoiding T, and we will see in the next 
section that these 3 variables can be combined in a more efficient way. The basic methods 
presented above become imprecise in some places of the kinematic domain. To have a single 
method covering the complete kinematic region measurable at HERA we can either make 
use of the kinematic fitting procedure (see sect. 5) or explore further combinations of the 4 
variables [E,9,T,,^] (sect. 3 and 7). 



3 The S method 

With 3 variables, we can determine y and independently of initial state QED radiation 
(ISR) by reconstructing the incident electron energy. From energy-momentum conservation, 
the measured quantity A is equal to two times the electron beam energy, if no particle escape 
detection: 

A = S + S(l - cos 0) = 2E'' (7) 

We can then obtain from y^ a new expression [§] which gives y at the hard interaction 
vertex even if an ISR photon is emitted. 

E + E[l — cosd) 

ys has another important characteristic: at high y, i.e in the most interesting region at low 
Q^, E dominates over E{1 — cos 6), and the fluctuations and errors on E start cancelling 
between numerator and denominator, allowing an optimal y reconstruction when using the 
hadrons. Writing the error propagation on yj], 

5y , , ,(5E 5E 59 , 
y h E tant^/2 

we can observe the following differences with y^ (eq.5) and yh {Syh/yh = (^E/E): 

i) At low y there is no divergence of 6y/y for y^,, which has practically the same behaviour as 
yh since the 5E/E term is experimentally bigger than the two terms coming from the electron 
measurement. 

ii) at high y, the error on y decreases as (1 — y) for but is still dominated experimentally 
by the error on E. The (^E/E term is partly compensated by the (1 — y) term contrarily to 
the yh case. 

In fig.l is shown the behaviour at low x-low (see sect. 4.1 for the description of the 
sample of Monte Carlo events used here) of the correlation / He compared to yh / ye and to 
UDA/Ue- These predicted behaviours have been reproduced by the HI collaboration for the 
F2 analysis on the 1993 data |Q. We can observe the expected improvement of at high y, 
despite the worsening of yh due to the low energy of the hadronic final state|^. In this plot 
y = 0.5 corresponds in average to x ~ 2.10"^. 



■^This effect is amplified in tfie Hf detector by tfie poor fiadronic calorimeter coverage in the backward 
direction, in which the hadronic final state particles tend to go at high y. We can remark that ys is relatively 
insensitive to this peculiarity. 
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Figure 1: Monte Carlo comparison of y reconstruction for several hadronic methods w.r.t. 
to Ue, at high (0.2-0.5) and very high y (0.5-0.8). Both bias and resolution improve with y 
for the S method although these become worse for the hadrons-only method. yoA has larger 
biases and resolutions than y^,, as can be read from the ht results written inside the hgure. 



Prom E, 9 and S we determine Q^, also independent of ISR when assuming that an initial 
state radiated photon which is colinear to the beam line do not carry transverse momentum. 

2 E"^ sin2 d 



Qi = ^— (10) 



From = xys we can obtain in two ways: 

i) taking the "true" s at the hard interaction vertex, as will be dicussed in section 7. ii) 
taking s as computed from the beam energies (s = AE^E^): 

Thus we have defined a new method, which has similar characteristics as the mixed method 
at low and medium y, but a higher precision at high y allowing to cross-check the e method 
in this region. It must be stressed that this improvement does not come from the precision 
of the electron measurement, on which on the contrary relies the precision of the e method. 
These methods then provide two measurements essentially independent. Finally y^, and Q^^ 
being ISR independent, the radiative corrections to the S method are smaller than to the 
other methods (see sect. 6). 



4 Comparison of Electron, DA, Mixed and E methods 



We study in this section the resolutions and biases of the x and variables as reconstructed 
by the different methods. The Monte Carlo event sample on which this study is based has been 
generated according to cross-section above 5 GeV^, with the DJANGO [10| program which 



relies on HERACLES |11] for the treatment of QED radiation, and for the simulation of the 



hadronic final state on the color dipole model (CDM) as implemented in ARIADNE [12|. The 
CDM is known to give the best description of the energy fiow on the HERA data The 
MRSH 1 15] parton densities parametrizations, which include the constraints from the 1992 
HERA data were used. The generated events were then processed through the HI analysis 
chain, which includes a detailed simulation based on GEANT |l^ of the HI detector |17]. 
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Thus the results presented here after include in a realistic way the effect of inhomogeneities 
or small miscalibrations encountered in a real detector. The conclusions that we will draw 
are nevertheless general, since the characteristics of the HI and ZEUS detectors are rather 
similar after the major detector effects are corrected by the reconstruction program. To 
make a meaningful comparison, all methods are applied on the same events using the same 



HIKINE [18 1 software package which is part of the HI reconstruction program. One of its 
characteristics is to use a combination of track momentum and calorimetric energy in the 
central part of the detector (25° < 9 < 155°) to improve the precision on the hadronic 
final state measurement and in particular to render it sensitive to the low energy hadrons 
{Efi < 0.5 GeV) present in large number at low Q^. This treatment is applied consistently 
to all reconstruction methods. 

The structure function analysis is generally done in bins of x and Q^, thus we will study 
here the following distributions^: x^ethodlx and Q'^ethodl Q'^ ■ h\as on the reconstructed 
variable is observed when the distribution peaks at 1. In order to conserve a good proportion 
of events originating from a given bin in that bin (defined as the "purity" of the bin) we need 
to have small biases and good resolutions on x and . The better the resolution the finer 
binning we can reach, assuming we have enough events, and thus the more subtle effects in 
the i*2 evolution along x or can be observed. This is particularly true at low x where the 
number of events collected in the first three years of operation is large enough to make the 
finest binning allowed by the detector resolutions. 

We will focus successively on two very different kinematic regions, for obvious physical 
reasons: the low x domain (x < 10^^) which at HERA is accessible only below ~ 50 GeV^, 
and the high domain which we will define as > 200 GeV^. In both cases, the events 
must also satisfy the following conditions: 6 < 173°, E > 8 GeV, A > 30 GeV, this last 
condition being helpful to reduce radiative effects and, on data, the background to the deep 
inelastic interactions. However for the distributions shown in the next two sections, only 
events without emitted photons from the electron line are used, in order not to mix the 
smearing due to the reconstruction method with the one coming from the radiation effect. 



4.1 Comparison at Low x 



For this study we ask: > 7 GeV^. The results are shown in figs. 2 and 3 for x and 
respectively, with one column per method and 5 bins in y. We distinguish very high (0.5- 
0.8), high (0.2-0.5), medium (0.1-0.2), low (0.05-0.1) and very low y (0.01-0.05), reminding 
that the lowest x are at highest y and viceversa. From these distributions we can draw the 
following conclusions which are valid in the low region: 



The S method has, as expected, a much more precise x reconstruction than the mixed 
method at high and very high y. The behaviour of xs and Xm is almost identical at 
low and medium y. 



*When X, y or have not the subscript which corresponds to a given method, they are assumed to be the 
true i.e. generated values at the hard interaction vertex. 
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Figure 2: Comparison Xme^/^od/^ ^ow (Q^ > 7 GeV^) for the E, mixed , DA and 
e methods. From top to bottom, each row represent a bin in y: very high (0.5-0.8), high 
(0.2-0.5), medium (0.1-0.2), low (0.05-0.1), very low (0.01-0.05). 
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Figure 3: Comparison Q'^methoJ Q"^ (Q"^ > GeV"^) for the S , DA and e methods. 

From top to bottom, each row represent a bin in y: very high (0.5-0.8), high (0.2-0.5), medium 
(0.1-0.2), low (0.05-0.1), very low (0.01-0.05). 
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• xj: is the only x which is reconstructed with small biases (< 10%) and precise resolutions 
(< 20%) over the full y range (0.01-0.8). xda can compete at low and medium y, but 
becomes very imprecise when y rises. Xe is slightly more precise than xs at high y, 
then rapidly becomes imprecise with decreasing y. 

• QI has always very small biases and optimal resolutions (~ 4%). has no bias 
and good resolution (~ 7%) at low and medium y, but becomes more imprecise with 
increasing y. Q'jjji has the same tendency, but in a more extreme way: better resolution 
at low y (~ 3%) and poorer resolution at very high y. 

In conclusion, for the low x physics , i.e. at low Q^, the E method allows to make a complete 
measurement in the full kinematic range. The e method remains the most precise method 
at high y, but now it can be compared to a hadronic method in the low x domain. The DA 
method is clearly disfavoured for low x physics at low Q^. Finally the mixed method can be 
advantageously replaced by the S method. 

4.2 Comparison at High 

For this study, we request > 200 GeV"^ and the Xmethod/x distributions are shown in fig. 4. 
The corresponding distributions have all a gaussian shape, so we just give their mean 
value and standard deviation in table 1. 









Qi/Q' 


0.5 < y < 0.8 
0.2 < y < 0.5 
O.K y < 0.2 
0.05 <y< 0.10 
0.01 <y< 0.05 


0.89 a = 0.16 
0.98 a = 0.08 
1.01 a = 0.06 

1.01 a = 0.05 

1.02 a = 0.05 


0.99 a = 0.07 
0.99 a = 0.03 
0.98 a = 0.03 
0.98 a = 0.02 
0.98 a = 0.02 


0.99 a = 0.04 
0.99 a = 0.03 

1.00 a = 0.03 
0.99 a = 0.03 

1.01 a = 0.03 



Table 1: Comparison Ql^ethod/Q^ ^^S^ (Q^ > 200 GeV'^) for the S , DA and e 
methods. 

Prom fig.4 and tab.l we can draw the following conclusions: 

• The general behaviour of the E, the mixed and the e method have not changed, both 
in X and Q^. A slight improvement in precision is visible at higher Q^, which is clearly 
understood by a more precise energy measurement of the electron and of the hadrons. 

• The DA method becomes more precise and improves at higher Q^: 

- In X it has better resolutions than the S method for y between 0.05 and 0.3. It is 
more affected by losses in the beam pipe at very low y, and not as precise as the e 
or S methods at high y, but this last feature will disappear at very high (above 
2000 GeV^). 

- In it is very precise at medium and low y, and display a similar behaviour as Q| 
at higher y, but without biases and better resolutions. The Q%y^ resolution remains 
below 10% over the full y range. 
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In conclusion, it is also possible to make a very precise kinematics reconstruction at high 
although this precision is not as important as at low due to the fast decrease of the DIS 
cross-section with Q^, which will always prevent at HERA a very fine "binning". 




soo 

400 

300 

zoo 

/ GO 
800 




petxA: s-igrrii 

0.99 o.oe ' 



peafc s-iffma 
0.99 O.tS 




jjea/c s'igma 
0.90 0.3S 




I ' I I I I I I > I I 



-pea-h s-ig-rrtcL 
O^G7 0^45 




Figure 4: Comparison Xmethod/x at high Q"^ (Q"^ > 200 GeV"^) for the S, mixed , DA and 
e methods. From top to bottom, each row represent a bin in y: very high (0.5-0.8), high 
(0.2-0.5), medium (0.1-0.2), low (0.05-0.1), very low (0.01-0.05). 
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4.3 Effect of QED Radiation on Kinematic Reconstruction 



In this section we study the QED radiation effect on the kinematic. We take the same 
selection criterias than in sect. 4.1, but requesting now the presence of a radiated photon 
having an energy such that 1.5 GeV < < 12.5 GeV, the upper bound coming from the 
A cut. The true x,y and are defined at the hard interaction vertex taking into account 
the ISR or final state radiation (FSR) effects: For the ISR, these variables are computed 
assuming an incident electron energy equal to — £">' . For the FSR, the photon is summed 
with the electron. This is experimentally happening in most of the events since the colinear 
photon is merged in the cluster of energy deposited by the electron in the calorimeter. 

Since the results are in general very similar to the ones obtained excluding the QED 
effects, we can just summarize this comparison on the biases and resolutions of the distribu- 
tions with and without radiation. 

i) the reconstruction is essentially not affected by radiation in all the methods studied. 
The additional biases and enlargement of resolution are always smaller than 2%. 

ii) the X reconstruction in the hadronic methods (m,DA,S) is also very stable against radia- 
tion, with changes of the order of 1 to 3%. 

iii) the most affected variable is Xg when going at low y. The additional bias due to radiation 
on Xe/x is in the 5 bins, from very high to very low y: —1%, —5%, —10%, —20%, —35%. 



5 Beyond the E Method 



5.1 The eE Method 



From the previous section, we can conclude that the S method can replace the mixed method 
and provide, with a single prescription, a precise F2 measurement over the complete range 
covered by at least one the 4 basic methods It is however clear from figs. 2-4 than the two 
most stable reconstruction for x and all over the kinematic plane are and Q^. We can 
then simply define a so-called eS method, in which the is determined from the electron 
quantities, and x is reconstructed according to the S prescription: 

2 E^ S 

xe^=x^ Ql^ = Ql y^^= i^ + Eil-cose)r ^''^ 

The improvement brought by this method can also be observed on y^Y;, by comparing its 
correlation ye's/ye to the ones shown in fig.l, as can be also read in table 2. The correlation 
with ye is excellent at high y, in fact higher than its correlations to the generated y. We 
thus have a method which optimizes the reconstruction quality, but has the disadvantage of 
becoming correlated to the e method at high y and also to enlarge the small influence of 
radiative correction displayed by the S method as we shall see in section 6. 



10 





0.2 <y < 0.5 


0.5 < y < 0.5 


Vh/Ve 
VDA/Ve 

ys/ye 
yeT./ye 


0.88 a = 0.26 
0.90 a = 0.22 
0.94 a = 0.17 
1.00 a = 0.06 


0.85 a = 0.27 
0.96 a = 0.16 
0.97 a = 0.11 
1.01 a = 0.03 


yT./y 
yenly 
ye/y 


0.95 a = 0.15 
0.99 a = 0.09 
0.99 a = 0.06 


0.95 a = 0.11 
0.99 a = 0.06 
0.99 a = 0.03 



Table 2: Comparison of y reconstruction quality, for the S, mixed, DA, e and eS method. 



5.2 Kinematic Fitting 



The ultimate improvement for the {x,Q'^) determination will come from a kinematical fitting 
procedure which use the full information in an optimal way. The weights on the redundant 
information coming from the electron and the hadronic final state have to be assigned de- 
pending on the kinematic of the event itself, and need, for a precise determination, a deep 
knowledge of the detector response. So far the fit has been performed using the conservation 
of A and of the transverse momentum [Q, |2^ (for simplicity, we do not consider here the 
staightforward modification needed to treat the ISR case). This is equivalent to say that the 
fit has been performed using [E, 9, S, T] as measured variables and in fact using only [E, 9, S] 
since T is badly measured. It is then a posteriori clear why the fit behaves only slightly better 
than the e method at high y, the S method at medium y, and the S or m method at low y. 

To go beyond this "3-variables fit" we need to include in the fitting procedure the infor- 
mation contained in the hadronic transverse momentum. This can be achieved by replacing 
(cf sect. 2) T by the hadronic angle 7. We can then fit y and according to the following 
equations 

E(l -cos 6*) = 2^*^(1 -y) Esm9 = jQ'^{l-y) (13) 



E = 2E=, tan2 = 2Ey^j^ (14) 

hence the constraint equations can be written as 

i: + E{l-cos9) = 2E^ tan2= ^^^-f.^7°^^) (15) 

2 E sni 9 

Assuming that the error covariance matrices has been determined for eqs. 13-14, this fit will 
then provide everywhere a better result than any of the e, h, m, DA, S methods, i.e. it will 
also have the properties of the DA method like its high precision at large Q^. 

A complete study of the kinematic fitting is beyond the scope of this paper but we 
can conclude that there are at least two ways to become more precise than the basic or S 
methods, which both have the drawback that the measurement they can provide cannot be 
cross-checked a posteriori in a simple and independent way. 
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Figure 5: Radiative corrections (excluding vacuum polarisation correction) as a function of 
X in two bins of for the S, eS, DA, e and mixed method, requesting two minimal scattered 
electron energy: E > 4 GeV (upper curve of each method) and E > 8 GeV (lower curve). 
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6 Comparison of QED Radiative Corrections 



The size Q'^) of the QED radiative correction (r.c.) which has to be apphed to the 

measured differential cross-section to obtain F2 depends on the method used, since the cross- 
section is determined in bins of reconstructed x and Q^. We have 

'''"'(2-2y + -^)F2(x,Q2) (16) 




where R = aL/aT is the ratio of the cross-section of longitudinally and transversely polarised 
photons. The r.c. due to the virtual graphs, which are dominated by the fermion loops in the 
self-energy of the photon (vacuum polarization) can be included in a "running" a = a{Q'^) 
in this equation |21]. These corrections are a function of only, ranging from 5 to 13 %, for 
Q2 between 5 and 10^ GeV^. The remaining Sj^c is dominated by the ISR from the lepton line 
and is thus a sensitive probe to distinguish the r.c. behaviour for the different reconstruction 
methods. 

The radiative corrections for the e, m,DA,S and eS methods were obtained with the 
HELIOS program[^] in which they are calculated in the Leading Log Approximation up to 
0{a'^), and cross-checked with the HERACLES program as implemented in DJANGO ||l^ 
for the e and S methods [Esl 



In fig.^ is shown as a function of 2; in 2 representative bins of (8.5 and 200 GeV"^) the 
size of the radiative corrections for every method, without cut on A. We nevertheless take 
into account the minimal electron energy required for DIS identification. Two typical values 
have been used, E = 4: GeV and E = 8 GeV, so two curves for each method are displayed. 
The variation of the r.c. for the 2 minimum energies can reach large values {S^/^ ~ 16%) for 
the DA or the e method, while it rather small for the m and eS ((J™'*^^ — 2%) and almost 
null for the S method. In the F2 measurable region (0.01 < y < 0.8) there is no dependence 
on X or neither for nor, at the 2% level, for Sf.f, contrarily to and Sf.^ (dependence 
in x) or (dependence in Q^). We can also remark that using the nominal s in the S 
method has a very small effect on the radiative corrections, since (5,^ ~ —1.5%. 

The radiative corrections are therefore experimentally very difficult to control at the 1% 
level for the e and DA methods without introducing a cut on A of the events e.g. A > 30 GeV 
which reduces the ISR effects. After such a cut all "hadronic" methods have r.c. displaying an 
almost constant behaviour in x and Q^, and amounting at 8.5 GeV'^ to 5*^* = drc + ^rc'^'^°^' — 
+2, 4-5, +6, -|-9 % for the S, eS, m, DA methods respectively. There is still a dependence on 
X for 5^^, which stops growing for y values above z = A/2E'^. 

We can conclude that the S method has small r.c. and that they barely depend on x, 
or the presence of the A cut, therefore allowing to make a structure function measurement 
free of radiative corrections uncertainties. It is again complementary to the electron method 
which has r.c. varying as a function of x, and also dependent on the parton densities in the 
proton in the absence of a A cut. 
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Figure 6: Comparison of x and reconstruction at low (Q'^ > 7 GeV^) for S and DA 
method independent of initial QED-bremsstrahlung (IE, ID A) From top to bottom, each row 
represent a bin in y: very high (0.5-0.8), high (0.2-0.5), medium (0.1-0.2), low (0.05-0.1), very- 
low (0.01-0.05). 
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7 DA and E Methods Independent of QED radiation 



In case of initial state radiation of a photon from the electron we have stme = ^{E"^ — E'^)E'p . 
To fully constrain the kinematics we need three independent variables and from the four 
variables defined in section 2, four combinations are possible. In practice only two ^ will 
deserve here our attention: the two which come directly from the DA and the S method. 

The first of these methods was developped in Q, by determining the initial energy E^ 
from [£',0,7]. Once E^ is reconstructed, it can be injected in the DA method to obtain the 
formulae for the double-angle method independent of ISR, labeled IDA here after: 

_ tan 7/2 + tan 61/2 _ E cot 7/2 + cot 61/2 

~ cot 0/2 + tan 0/2 ^^^^ ~ E^ cot 0/2 + tan 0/2 

^9 ^9 tan7/2 + tan 0/2 

Q]oA = ta„«/2 (18) 



The other method is the S independent of ISR (IS), which is obtained by just using str 
instead of s [24|. We have y/s = ys, Q/s = ^^'^ 



Ej-B 



J: + E{1-cos9) 



E cos^e/l 
E^ yn 



(19) 



A comparison of the IS and IDA methods is presented in fig. 6, following the scheme of figs. 2-3. 
The following conclusions can be drawn. 



• At low Q^, the IS method is slightly more precise than the IDA one in all y bins both 
for X and for Q^. 

• xiY, is less precise than xy. at high while they have similar resolutions at low and 
medium y. 

• xiDA and Q^dj^ are more precise than their DA counterpart at low x, but have wider 
resolutions at low y. 



The loss of precision when trying to reconstruct E^ is compensated by the reconstruction of 
the true kinematics when studying the special class of events with catastrophic initial state 
radiation. However for the structure function measurement it is better, as it has been done on 
the 1992 and 1993 data by the HI and ZEUS collaborations 0, |, ||, ||], to use the methods 



described in the preceding sections (e, m, DA, S) since the additional smearing observed in 
these two ISR independent methods is larger than the smearing introduced by QED radiation. 
These ISR independent methods can also be used as has been demonstrated by a complete 
F2 analysis based on the IS method which has been done inside the HI collaboration on the 
1992 data |27], in good agreement with the published measurements. 



^The other two methods uses 2 variables from the hadronic final state and only one from the electron, thus 
being less precise. 



15 



8 Conclusion 



We presented a new prescription for DIS kinematic reconstruction well adapted to study the 
low x-low physics in the full y range (0.01 < y < 0.8) and which allows to measure F2 in the 
complete kinematic domain accessible at HERA. This S method can be modified to become 
fully independent of ISR (IS method) although this is not necessary, since the radiative 
corrections for the E were shown to be already very small, and some loss of precision when 
using the IS method was observed. By replacing by we defined another method 
which further optimizes the reconstruction precision at the cost of becoming correlated to 
the e method at high y. For the kinematic fitting we recommended to consider the hadronic 
angle instead of its transverse momentum as the second variable (beside S) characterizing the 
hadronic final state. After this review on the subject of reconstruction of the deep inelastic 
scattering kinematics at HERA, we expect that future improvements in this field will mostly 
come from deeper knowledge and better use of the detector response. 
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Appendix: 



y, and x formulae 



method 



X 



m 



DA 



IDA 



IS 



S 
2E^ 



Vh 

tan 7/2 
tan7/2+tan6»/2 

E+E(l-cos6i) 
VDA 



4^e2 



cot 61/2 



tan 7/2+tan 61/2 
sin2 ^ 



2 , ^„ e tan7/2+tang/2 
^^^2cot6'/2+tan6l/2 



QVys 
QVys 
QVys 
QVys 

QVys 

E cot 7/2+cot 61/2 
EP cot 6'/2+tan 61/2 

E_ cos^ g/2 
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